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ABSTRACT 

Observations of the cyclotron resonance scattering feature in the X-ray spectrum 
of GX 301-2 suggest that the surface field of the neutron star is -Bcrsf ~ 4 x 10 12 G. 
The same value has been derived in modelling the rapid spin-up episodes in terms of 
the Keplerian disk accretion scenario. However, the spin-down rate observed during 
the spin-down trends significantly exceeds the value expected in currently used spin- 
evolution scenarios. This indicates that either the surface field of the star exceeds 
50-Bcrsf, or a currently used accretion scenario is incomplete. We show that the above 
discrepancy can be avoided if the accreting material is magnetized. The magnetic 
pressure in the accretion flow increases more rapidly than its ram pressure and, under 
certain conditions, significantly affects the accretion picture. The spin-down torque 
applied to the neutron star in this case is larger than that evaluated within a non- 
magnetized accretion scenario. We find that the observed spin evolution of the pulsar 
can be explained in terms of the magnetically controlled accretion flow scenario provided 
the surface field of the neutron star is ~ -Bcrsf- 

Subject headings: X-rays: individual (GX301-2/4U 1223-62/Wray 977) — stars: neu- 
tron — stars: magnetic fields — X-rays: binaries: close 



1. Introduction 



A strong magnetizatio n of the n e utron star in the High Mass X-ray Binary (HMXB) GX 301-2 
has first been suspected bvlLipunov (| 19821 ). Assuming the star to rotate at the equilibrium period 
given bv lDavidson OstrikerT (|l973h he limited its dipole magnetic moment to > 10 32 Gem 3 . This 
implies the surface field of the star to be in excess of 2 x 10 14 G. 
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This estimate has been ch allenged by observations of the cyclotron resonance scattering feature 
(CRSF) at ~ 35 — 45keV (see lLa Barbera et alj|2005l . and references therein). The field strength 
of the pulsar inferred from th ese observation s, -Bcrsf ~ 4 x 10 12 G, is a factor of 50 smaller than 
the surface field evaluated by iLipunovl ([1982). 



However, a questio n about the field strength of GX 301-2 has recently been raised again by 
Doroshenko et al.l (|2010i ) who reported the pulsar to spin down on a time scale of years at the average 



rate of P s ~ 4.25 x 10 8 ss 1 (z> sh = Pn/P* ~ —10 13 Hz s l ). Using the spin-down scenarios of 



Davidson Qstrikerl ()1973l ) and llllarionov fe Kompaneeta (|1990l ) they put the same lower bound 
on the surface field of the neutron star as ILipunovl (|1982l ). To avoid the apparent discrepancy 
they have suggested that CRSF is generated at the top of a radiative-dominated accretion column, 
which is extended above the stellar surface up to an altitude of a few stellar radii. 

In this paper we argue that this conclusion is premature. The spin behavior of the pulsar is 
very complicated. The long-term spin-down trends are superposed with phases of stable rotation 
and rapid spin-up episodes (see Figure[IJ and Sect.[2|). The value of torques applied to the star 
strongly depends on the geometry and physical conditions in the accreting material (Sect. [3]). The 
pulsar behavior during the rapi d spin-up episodes can be explained in terms of accretion from a 
Keplerian disk (|Koh et al.lll9971 ) provided the surface field of the star is ~ -Bcrsf (Sect.HJ). The 
spin-down rate of the star during the spin-down trends strongly depends on the conditions in the 
material surrounding the pulsar magnetosphere (Sect. [5]). We find that the neutron star can brake 
harder if the accreting material is magnetized. The accretion flow in this case is approaching the 
magnetosphere of the neutron star in a form of a dense magnetic slab confined by the magnetic 
field of the flow itself. The interaction between the slab and stellar magnetosphere leads to an 
effective angular momentum transfer from the star to the slab (Sect.[6|). We show that the observed 
spin-down rate can be explained within this magnetically controlled accretion scenario provided the 
surface field of the neutron star is about -Bcrsf (Sect. [7]). Assumptions adopted in this accretion 
picture are briefly discussed in Sect. [8] and our conclusions are summarized in Sect.[9l 



2. Basic parameters 

GX 301-2 is a 685 s accreting X-r ay pulsar in a 4 1.5 d eccentric orbit around the early B- 
type supergiant companion Wray977 (ISato et al.l Il986l ). The massive companion underfills its 
Roche lobe and looses material at th e rate M ni it ^ 10~ 5 M(^yr~ 1 in the form of a relatively slow 
(v w ~ 300 - 400 km s" 1 ) stellar wind feaoer et alJbood ). 

As the neutron star moves through the wind of density with a relative velocity v re i = 
v ns + v w it captures material at a rate 9JT < 9H C , where 



(1) 



Here R ns is the radius, M ns the mass and v ns ~ 250 km s is the orbital velocity of the neutron 
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star, and Rg = 2GM ns /v^ cl is its Bondi radius. 



The interaction between the accreting material and the stellar magnetic field leads to formation 
of a magnetosphere which, in the first approximation, prevents the accretion flow from reaching 
the stellar surface. The flow is decelerated by the stellar magnetic field at a distance r m , which is 
defined by equating the pressure of the accretion flow with the magnetic pressure due to the dipole 
field of the neutron star and is referred to as the radius of the magnetosphere. The accretion flow 
enters the pulsar field at the magnetospheric boundary and flowing along the field lines reaches the 
surface of the neutron star at the magnetic pole regions. The rate at which the accreting material 
enters the pulsar field, 9Jl; n , and, correspondingly, the mass accretion rate onto the stellar surface, 
S0T a , can be evaluated as 



2 x 10 37 ergs 



-i 



gs 



(2) 



where Rq = i? ns /10 6 cm, m = M ns /1.4MQ and Lx is the X-ray luminosity of the s ystem, which is 
normalized to its average value Lx ~ (1— 3) x 10 37 erg s _1 for a distance of 1.8-3 kpc, ( Chichkov et al 
19951 : iKaper et alJbood ). 



The persistent character of the pulsar indicates that the magnetospheric radius of the neutron 
star is smaller than its corotation radius, 



GM n 



jj- 



1/3 



1.3 



10 10 m 



1/3 p 2/3 
^685 



cm, 



(3) 



and hence, the centrifugal barrier at the magnetospher ic boundary does not prevent the accretion 
flow from reaching the stellar surface (IShvartsmanlll970l ). Here P685 is the spin period of the neutron 
star in units of 685 s, and u s = 2ir/P s is its angular velocity. 

The pulsar shows a complicated picture of spin variation (see Figure[T]), which can be roughly 
divided into three p hases. The spin-down t rends with the average rate |z>(jo| ~ 10 _13 Hzs _1 and 
1-3 years duration (jDoroshenko et al.ll2010l ) are supe rposed with the short-term (about 30 days) 



rapid spin-up episodes with z> u o ~ 5 x 10 Hzs (|Koh et al.l 11997), and the phase of almost 
stable rotation (y < 10~ 13 Hzs -1 ), which lasts up to a few years ( Bildsten et al.l Il997l ). The 
pulsar transition between these phases occurs without any significant changes in the intensity and 
spectral parameters of the X-ray emission. In particular, X-ray luminosity of the system during 
rapid spin-up epi sodes is only a factor of 2-4 larger than that observed during the spin-down trends 
(|Koh et al.lll997l ). This indicates that the mass-transfer between the system components operates 
almost stationary (i.e. S9T C ~ 9?T a = SOT) without any significant changes in the parameters of the 
accretion flow. 
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3. Accretion-driven spin evolution 

The equation governing the spin evolution of an accreting neutron star reads 

Iu B = K su - K sd . (4) 

Here I is the moment of inertia of the neutron star, and K su and K sd are the spin-up and spin-down 
torques applied to the star from the accretion flow. 



3.1. Spin-up torque 



Since the material captured by the neutron star in a binary system possesses specific an- 
gular momentum, the mass accretion onto the neutron star is acco mpanied by the accretion of 
angular momentum at the rate J = £ f2 or b Rq Tl c (for discussion see, llllarionov fc Sunvaevl Il975l ; 
Illarionov fc Kompaneetslll99dl . and references therein). Here ^ rb = 27r/P or b is the orbital angular 
velocity and £ is a parameter accounting for dissipation of angular momentum in the accreting 
material (see e.g. lRuffertlll999i . and references therein). 



The geometry of the accretion flow in this case deviates from spherical symmetry. The devi- 
ation, however, will be significant only if r c i rc > r m , where r c ; rc = J 2 /GM ns 9Jt^ is a circularization 
radius, at which the angular velocity of the accr eting material, aw. = 6 ffi orh (Rq/'t) 2 , reaches the 



1 /9 

Keplerian angular velocity, lo^ = (r 3 /2GAf ns ) ( Bisnovatyi-Koganl Il99ll ) . If this condition is sat 



isfied a formation of the Keplerian disk can be expected. Otherwise, the accretion should be treated 
in the Quasi-Spherical (QS) approximation, which implies < co ea (r m ) <C Wk(r m ). 

The radius of the magnetosphere of a neutron star accreting from a non- magnetized wind 
(validity of this approximation is discussed in Sect. [6]) is ?" m ~ r a , where (see, e.g. lAronslll993l . and 
references therein) 



\VJl(2GM n 



a/2 



2/7 



> 



(5) 



r a for v re \ one 



Here \i is the dipole magnetic moment of the neutron star. Solving inequality r c 
finds that the neutron star in a wind-fed HMX B would be surrounded by a persistent Keplerian 
disk if v Te i < v cr , where (see e.g., Ilkhsanovll2007l ). 

-1/4 

(6) 



200 £&30 1/14 ™ 



ll/28oj>l/28 / -forb 
' 17 1 41.5 d 



kms 



Here ^30 = /i/10 30 Gcm 3 , 5^17 = 9Jt/10 17 gs 1 and £0.2 = £/0.2 is normalized to its average value 
deriv ed in numeric al modeling of quasi-spherical accretion in the non-magnetized flow approxima- 
tion (jRuffertl Il99gl . and references therein). The specific angular momentum of the material in 



the Keplerian disk is jk( r ) = UJ 'kT 2 - Therefore, the spin-up torque associated with the accretion of 



mate rial from the inner radius of the disk is Ksu = 9^ik(^m)^m = {GM ns r m ) 1/z (jPringle fc Rees 
19721 ). 
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Thus, the spin-up torque applied to an accreting neutron star in a wind-fed HMXB can be 
expressed as 

' Tl (GM ns r m ) 1/2 for v rel < v CT (disk) 

^su * I (7) 

971 £ il orb Rq for t; re i > v CT (QS) 



3.2. Spin-down torque 

The spin-down torque applied to an accreting neutron star is associated with interaction be- 
tween its magnetic field and material surrounding the magnetosphere. This interaction leads to a 
distortion of the magnetospheric field, turbulization of the material at the magnetospheric bound- 
ary, excitation of MHD waves and magnetic reconnection. The efficiency of these processes depends 
on the field geometry, conductivity of the accreting material and spectrum of the turbulent motions. 
A detailed modeling of these processes is, therefore, very complicated and is beyond the scope of this 
paper. Here we focus on a simplified task targeting at evaluation of upper limits to the spin-down 
torque under various assumptions about the geometry and parameters of the accretion flow. 

The spin-down torque applied to a neutron s tar a ccreting material from a Keplerian disk 
has first been evaluated by iLvnden-Bell Pringld (| 19741 ). The stellar magnetosphere in this case 
contains both closed and open field lines, which are extended beyond the inner radius of the disk. 
The angular velocity of the material in those parts of the Keplerian disk, which are located inside 
the corotation radius, exceeds the angular velocity of the star itself. The interaction between these 
parts of the disk and stellar magnetosphere tends, therefore, to spin-up the star. The angular 
velocity of the material in the Keplerian disk is, however, decreasing with radius and becomes 
smaller than the angular velocity of the neutron star in the region r > r cor . The interaction 
between the magnetospheric field lines and the disk at this distance tends to spin-down the star. 
As shown by ILvnden-Bell fe Pringld (|1974l ). the spin-down torque associated with this interaction 



can be evaluated as K$ = kt/j? /r^ or , where kt < 1 is the efficiency parameter. 

The magnetosphere of a neutron star undergoing spherical accretion is closed. The stellar 
field is screened by the surface currents at the magnetospheric boundary and there is no field 
lines extended beyond the magnetospheric radius. The boundary is convex towards the accreting 
material and is closed to the cusp poin t s located at the stellar magnetic axis (for discussion see 
Arons &: Lealll976l ; Eisner Sz Lamblll977l ; lMichellll977l ). This indicates that the interaction between 
the stellar magnetosphere and the accretion flow occurs solely at the the magnetospheric boundary. 
For the neutron star in this case to spin-down the angular velocity of the accreting material at the 
magnetospheric boundar y, w en (r m ), should be smaller than the angular velocity of the star itself, u s 



(jBisnovatyi-Koganl Il99ll ) . We assume here that this condition is satisfied and discuss the validity 



of this assumption in SectEJ 



The spin-down torque applied to a neutron star accreting form a quasi-spherical accretion flow, 
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in the first approximation, can be evaluated by finding the torq ue on a rotating sphere of the radius 
r m in a viscous medium (for discussion see e.g.. lLipunovlll992l ). The spin-down torque in this case 
can be expressed as 

K sd = ATirl^tpQV^, (8) 

where V\ is the viscosity coefficient, po is the density of material at the magnetospheric boundary 
and = [uj s — uj en (r m )} r m is the azimuthal component of the relative linear velocity between the 
sphere and the surrounding material. For a simplicity we consider the case of w en (r m ) <C uj s and 



hence, 
torque. 



uj s r m . This approximation is reasonable for evaluation of upper limits to the spin-down 



Assuming a turbulent nature of the viscosity, u t = hvtit, and taking into account that the 
density of the free-falling material at the magnetospheric boundary of the neutron star is po = 
9Jl c /47rr^ 1 Vff(r' m ), one finds 



K sd = kt lo s r m £ t (r m ) Wl c 



vs(r„ 



0) 



Here v t (r m ) and £t(r m ) are the velocity and the scale of turbulent motions, vg(r m ) = (2GM ns /r m ) 1 ^ 2 
is the free-fall velocity at the magnetospheric boundary and kt < 1 is the efficiency parameter. If 
the turbulent motions in the material at the magnetospheric boundary are excited by its interaction 
with the rotating magnetosphere, the viscosity coefficient can be expressed as v^' = ktVj,r m . This 
implies that the scale and the velocity of turbulent motions are limited to it < r m and Vt < va,. 
Putting these parameters to Eq. ([9]) and setting r m = r a one finds 

,2 



K. 



(o) 

sd 



V 2 

GMns 



k t — . 

'"cor 



(10) 



This upper limit t o the spin-down torque has previously been evaluated by iLipunovl (|1982l ) and 
Bisnovatvi-Koganl ()199ll ) for the case of quasi-spherical accretion. 



The neutron star under the same conditions is braking harder if u t > . This can be realized 
if Vt > (note, that the scale parameter, it, has already been limited above to its maximum 
possible value). On the other hand, the velocity of turbulent motions is limited to the speed 
of sound in the accreting material since the supersonic turbule nce it effectively suppresses by the 
Landau damping (for discussion see, e.g., Davies Pringle1ll981 ). Therefore, the maximum possible 
value of f t can be found by setting u t ~ vg, i.e. assuming that the velocity of the turbulent motions 
in the accretion flow is close to the speed of sound in the material heated in the shock at the 
magnetospheric b oundary up to the adiabatic temperature (for discussion see, e.g. lArons fc Lea 



1976 



Eisner & Lamb! 1 19771 ). Putting this to Eq. @ yields 



K. 



(t) 



fet Tl c w s r B 



L sd — "t -" l c^s ' m- (11) 

Finally, taking into account that = p? / (2r^GM ns )^ 2 (see Eq.[5|), one can express the maximum 
possible spin-down torque as 



K. 



(t) 

sd 



l£_ 

(r m r CO r) 



3/2 " 



(12) 



-7- 



The condition v t > can be satisfied if the neutron star accretes m aterial from eithe r a ho t 
turbulent atmosphere (as it occurs in the subsonic propeller scenario of iDavies &: Pringlelll98l] ). 
or a magnetized flow with a high value of the magnetic viscosity. The subsonic propeller scenario 
can be used, however, for interpretation of a limi ted number of X-ray pulsars. The spin period of 
subsonic propellers is limited to P s < P^ r , where (|Ikhsanovil200lh 



16/21 ffo-5/7 



P br ~ 450 ^ m 



15 



m 



-4/21 



(13) 



is a so-called break period. If the spin period of a neutron star exceeds Pbr the cooling of the 
atmosphere due to bremsstrahlung emission and turbulent motions becomes more effective than 
heating (due to propeller acti on by the neutron star). The X-ray source in this case would be 
observed as an X-ray burster (|Lamb et alJ 119771 ). Furthermore, the X-ray luminosity of subsonic 
propellers is limited to Lx < ^brj where 



Lb r ~ 10 34 H% 7 m~ 



- 2 / 7 2H- 4/7 



500 s 



erg s 



(14) 



Otherwise, cooling of the material at the magnetospheric boundary due to inverse Compton scat- 
tering of the X-ray photons, emitted from the surface of the neutron star, on the hot atmospheric 
electrons dominates heating due to propeller action. This condition has been derived by solving 
inequality t c (r m ) > t^, where 

3it r 2 m e c 2 



tc(r) 



propeller action (IDavies &; Pringle 



1981 



2a T L x 



(15) 



is the Compton cooling time (lElsner Lamb! 119771 ) and ~ l/w s is the heating time due to 



J). Here m e is the electron mass and <tt is the Thomson 
cross-section. Thus, the subsonic propellers can only appear as relatively low luminous, fast rotating 
pulsars, which is obviously not the case of GX 301-2. 

A high viscosity can also be expected if the material surrounding the neutron star is magnetized. 
The magnetic viscosity coefficient in this case can be expressed as v m = k t V\H. m , where Va = 
Bf/^irpo) 1 / 2 is the Alfven velocity and £ m is the characteristic scale of the magnetic field, Bf, in 
the material at the magnetospheric boundary. A rapid braking of the neutron star in this case can 
be expected if £ m ~ r m and Va > V^. The Magnetically Controlled Accretion scenario in which 
this situation is realized will be discussed in Sect.El 

Concluding this section one can express the spin-down torque applied to a neutron star from 
the accretion flow in the following form 



h — 

'"cor 



for f t < 



(0) 



(16) 



for i>t > f t 



(0) 



( \ 3 / 2 

y m ^*cor ) 

where kt < 1 is a parameter accounting for the conductivity, spectrum of turbulence and inhomo- 
geneities of the accreting material. 
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4. Rapid spin-up episodes 

The spin-up torque applied to the neutron star during the rapid spin-up episodes is K su = 
K + K sd , where 

Kq = 2ttIu su ~ 3 x 10 34 ^(^su/^uo) dyne cm, (17) 

and ^5 = //T0 45 gcm 2 . If the surface field of the neutron star does not considerably exceed 
-BcRSF the spin-down torque is significantly smaller than Kq. In particular, the spin-down torque 
applied to the neutron star from the Keplerian disk under these conditions (B* = £>crsf) is 
= /c t /" 2 / r cor ~ 2 x 10 30 /c t dyne cm (see Eq. [To]) . The spin-up torque in this case can be 
evaluated as K su > Kq. Solving this inequality with spin- up torque expressed by Eq. ([7]) one finds 
that the spin behavior of the pulsar during the rapid spin-up episodes can be explained in terms 
of accretion from the Keplerian disk provided the dipole magnetic moment of the neutron star is 
M > Mo, where 

/ ■ \ 7 1 2 

/ u ^2xl0 30 Gcm 3 / 4 7 5 /2 93tr7 3 ^ 3/2 ( — ) • (18) 
This corresponds to the surface field Bq = 2/io/-R 3 s >4 x 10 12 G, which is close to I?crsf- 

An attempt to explain the spin-up episodes in terms of the QS accretion scenario encounters 
difficulties. The spin-up torque in this case could be as high as Kq only if the relative velocity of 
the neutron star were v re \ < vq, where 



200 ^ - 1/2 < 4 (j^y^kmB- 1 . (19) 



v 41.5d 



However, a formation of the Keplerian disk under these conditions cannot be avoided (see Eq. [6J. 
It therefore appears that the interpretation of the spin-up episodes in terms of the Keplerian disk 
accretion scenario is reliable and consistent with the surface field of the neutron star measured 



through observations of CRSF (for discussion see also iKoh et al.1 119971 ) . 



5. Spin-down trends 

The spin-down torque applied to the neutron star during the spin-down trends is \K sd \ = 
\Ki\ + |-f£T S u|, where 

LKi| = 2vr/z> sd ~ 6 x 10 32 J 45 ( — ] dyne cm. (20) 

An interpretation of the pulsar behavior during these trends cannot be made within the existing 
spin-down scenarios without additional assumptions. If the star accretes material from a Keplerian 
disk the spin-up torque is |-fC s ^P| ^> \K\\ (see above). The spin-down behavior of the pulsar in this 
case implies \Kj\ '\ > \Kbu |- To satisfy this condition one has to assume that the surface field of 
the star is in excess of 5 x 10 G, which is 100 times larger than i?cRSF- 
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The spin-up torque applied to the neutron star is significantly smaller if the star undergoes 
quasi-spherical accretion. Its value depends on the efficiency of the processes responsible for angular 
momentum dissipation in the accreting material which is accounted for by the parameter £ (see 
Eq. [7]). This parameter in the case under consideration is limited by the condition cj en (r m ) < u s , 
which implies that the angular velocity of the accreting material at the magnetospheric boundary is 
smaller than the angular velocity of the star itself. Otherwise, the interacti on between the accretio n 



flow and stellar magnetosphere tends to spin- up the star (for discussion see. iBisnovatyi-Kogan 
Setting r m = r a and solving the above inequality for the parameters of GX 301-2 yields 



1991 



£ < 



0.03 m-W L- 7 ^ Rf /7 (^) x 

-Y(—) 



(21) 



8/7 



400 km s 



This indicates that | during the spin-down trends does not exceed \K\ | and hence, the condition 
for the spin-down phase can be expressed as |K s d| > \K\\. 

Putting K s d ~ Kjty one finds that the above condition can be satisfied only if the surface field 



of the neutron star is > 10 14 k t G (for discussion see also lLipunovlll982l ; lDoroshenko et al.ll2010l ) 



To satisfy the condition > K\ one has to assume that either the surface field of the neutron 



— 7 /X 

star is > 3 k t -Bcrsf , or the radius of the stellar magnetosphere is smaller than its canonical 
value, r a , at least by a factor of 2k t 



-2/3 



Thus, both the Keplerian disk and quasi-spherical accretion scenarios encounters difficulties ex- 
plaining the spin behavior of GX 301-2 during the spin-down trends. This may indicate that either 
the surface field of the star is indeed in excess of -E>crsf> or the neutron star in this binary system 



accret es in a different way. The first possibility has been already discussed by iDoroshenko et al 



(|2010l ). Here we focus on the analysis of an alternative accretion scenario. We show that the 
peculiar spin behavior of GX 301-2 can be explained in terms of the magnetic accretion scenario 
(|Shvartsmanlll97ll : IBisnovatyi-Kogan Ruzmaikinlll974l . Il976l ) provided the material captured by 
the neutron star from the wind of its normal companion is magnetized. 



6. Magnetic accretion 



Let us consider a situation in which the magnetic energy density in the material captured by 
the neutron star, £ m (i2c) = Bf /8iT, is comparable to its thermal energy density, £ t h(-Rc) = Ax>c 2 , 
i.e. (3 = £th/£m ~ 1- Here c s o is the speed of sound and Bfo is the strength of the magnetic field 
in the captured material at Rq. The magnetic field in the free- falling material is dominated by its 
radial component, B v , ( the transverse scales in the free- falling flow contract as r -2 , while the radial 



scales expand as r 1 / 2 . IZeldovich fc Shakur3ll96flh. which under the magnetic fl ux conservation 



condition increases as B T (r) ~ B i0 (R G /r) z (jBisnovatvi-Kogan fc Fridmajll970l ). The magnetic 
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pressure in the free-falling material, 



£m(r) = &m(R G ) 



Rg 



increases, therefore, more rapidly than its ram pressure, 

£ram(r) = £ ram (i? G ) (— j 



5/2 



(22) 



(23) 



and hence, the gravitational energy of the star is converted predominantly into the magnetic pres- 



sure of the free-falling material, £ m /£ ram oc r 
the captured material at the Bondi radius. 



- 3 / 2 . Here E iam {R 



G 



P™ V re\ 



is the ram pressure of 



The distance i? s h (hereafter Shvartsman radius) at which the magnetic pressure in the accretion 
flow reaches its ram pressure can be derived by solving equation £ m (i? s h) = £ram(-R s h)- This yields 
(|Shvartsmanlll97ll ). 

\ 4/3 



R 



ah 



-2/3 



VrelJ 



Rg = P 



-2/3 



2GM ns ct /3 

10/3 
V rel 



(24) 



If the magnetic flux in the accreting material is conserved, the accretion ends at the Shvartsman 
radius. Further accretion in this case is impossible. Otherwise, the magnetic energy in the flow 
wo uld exceed its grav itational energy, which contradicts the energy conservation law (for discussion 



sec 



Shvartsmanlll97ll ) . Therefore, the accretion flow can approach the star to a closer distance only 



if dissipation of the magnetic field in the flow occurs. If the field dissipation is governed by magnetic 
reconnection the characteristic time of the accretion process inside i? s h is limited to t > t rec , where 



VmVA 



"i 1 iff 



VA 



(25) 



The value of the efficiency parameter ry m depends on physical conditions and field co nfiguration in 



the region of reconnection, and in the general cas e is limited to < r] m < 0.1 (see, e.g. JParker 



1971 



Kadomtsevl[l987l : iNoglik et allhood : ISomovll^Ood . and references therein) . The Alfven velocity in the 
accreting material increases as it approaches the neutron star from the initial value VXo ~ /? _1//2 c s 
to VX(-Rsh) ~ Vff. Hence, the time of the magnetic reconnection in the accretion flow remains 
significantly larger than the dynamical (free-fall) time, % = r/Vg = (r 3 /2GM ns ) 1 ^ 2 , at any stage of 
the accretion process. This proves reliability of the assumption about conservation of the magnetic 
flux in the free-falling material. But, on the other hand, it also suggests that the accretion flow 
will be decelerated at the Shvartsman radius by its own magnetic field and the accretion process 
in the region r < R s h operates in the diffusion approximation. 

Rapid amplification of the magnetic field in the spherical flow and deceleration of the flow 
by its own magnetic field at the Shvartsm an radius have been confirmed in analytical studies 



( Bisnovatvi-Kogan &: Ruzma ikin 



1974 



1976) and num e rical calculations of magnetized spherical 



accretion onto a black hole (jlgumenshchev et al.l 120031 ; llgumenshchevi 120061 ) . These calculations 
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have shown that the magnetized flow is shock-heated at the Shvartsman radius up to the adiabatic 
temperature. The structure of the accretion flow inside the Shvartsman radius depends on the 
efficiency of cooling of the accreting material. If cooling is inefficient the accretion process switches 
i nto the convective-domina ted stage in which some material is leaving the system in a form of jets 
(jlgumenshchev et al.ll2003l ). Otherwise, the material tends to flow along the lines of the large scale 
field of the accreting matter itself an d accumulates in a dense non-Keplerian slab (see Fig. 1 in 



Bisnovatyi-Kogan &; Ruzmaikinlll976l ). The material in the slab is confined by its own magnetic 



field and its radial motion continues as the field is annihilating. This indicates that the accretion 
process in the slab occurs on the reconnection timescale expressed by Eq. 



6.1. Magnetic accretion in X-ray pulsars 

The Shvartsman radius exceeds the canonical magnetospheric radius of the neutron star evalu- 
ated in the quasi-spherical accretion scenario, r a , if the relative velocity of the neutron star through 
the wind of its massive companion satisfies the condition v re i < j; mca , where 

w m ca = /9~ 1/S (2GM ns ) 12 / 35 /i- 6 / 35 SOT 3 / 35 c s 2/5 . (26) 

For typical parameters of HMXBs this velocity is 

, mca ~ 680 /T 1/5 ™ 12/35 /%? /35 4 /5 kms-\ 

which under the conditions of interest substantially exceeds the critical velocity at which the forma- 
tion of the Keplerian disk in the system can be expected, v cr (see Eq. [6]). Here cq = c s o/10 6 cms -1 . 
This finding allows us to distinguish a subclass of HMXBs in which the accretion occurs in a spher- 
ically symmetrical fashion and the accreting material is strongly affected by the magnetic field of 
the flow itself. This subclass is defined by the condition v cr < v re \ < w m ca- 

The structure of the accretion flow inside the Shvart sman radius depend s on the efficiency of 



coolin g processes in the accreting material. As shown by lArons Leal (|1976l ) and lElsner &: Lamb 



(|1977l ). the cooling of the material accreting onto a neutron star is dominated by the inverse 
Compton scattering of X-ray photos emitted from the stellar surface on the hot electrons of the 
accretion flow. This mechanism in the considered case will be effective if the Compton cooling time 
of the accretion flow at the Shvartsman radius, i c (i? s h), is smaller than the characteristic time of 
the accretion process, which in the case of the magnetic accretion scenario is the reconnection time, 
t rec - Combining Eqs. (fT5j) and (i25j) . one can express inequality i c (-Rsh) < ^rec as Lx > L CT , where 

L cr * 3 x 10- ™' /2 ^ 1/8 (gfc) (^y 2 ^- 1 - <»> 

This indicates that the cooling of magnetized flow can be effective even in faint X-ray pulsars in 
which the Shvartsman radius does not significantly exceed the Alfven radius of the neutron star. 
The structure of the accretion flow in those systems in which the above conditions are satisfied can 
be treated in terms of the magnetic slab. 
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6.2. Magnetospheric radius 



The material in the slab is approaching the neutron star up to a distance at which its gas 
pressure becomes equal to the magnetic pressure due to the stellar dipole magnetic field. The gas 
density in the slab at this distance (the magnetospheric radius) is, therefore, 



Psl 



(j, 2 m p 



27r k B T 

where Tq is the temperature of the material at the inner radius of the slab. 



(28) 



The value of the magnetospheric radius in the general case depends on the mode by which 
the m aterial enters the stellar field at the magnetospheric boundary. As shown bv lElsner &: Lamb 
(j 19841 ). the diffusion rate of the accreting material into the stellar field can be evaluated as 



SKin(r-m) = 4:irr m 5 m p V s (r B 



(29) 



47rrU A D l e ( i 2 p (2GM ns ) 1 / A , 



where S m = (D c g r^) 1 ^ 2 is the thickness of the diffusion layer at the magnetospheric boundary 
(magnetopause) and D e Q is the effective diffusion coefficient. The diffusion time is determined 
by the time on which the material being penetrated into the field leaves the magnetopause by 
free-falling along the magnetospheric field lines towards the stellar surface, ~ ts(r m ). 



The estimate (|29l) has good theoretical and observational grounds. It has been justified by 
studies of the Earth's magnetosphere which show that the rate of solar wind penetration into the 
magnetic field of the Earth can be evaluated taking D c s ~ Db, where 



Db = CtB 



ck B T rl 
2e/u 



(30) 



is the Bohm diffusio n coefficient, e is th e electron charge and «b is the efficiency parameter, which 



ranges in 0.1 — 0.25 ([Gosling et al.lll99ll ). The diffusion process i n this case is governed by magnetic 
reconnection and drift-dissipative instabilities (|Paschmannll2008l . and references therein). Note also 
that the same diffusion rate has previously been measured in plasma experiments with TOKAM AKs 



(jKadomtsev &: Shafranovlll983l ) and evaluated from observations of solar flares (|Priestlll982l ). 



A higher rate of plasma penetration into the magnetosphere could be expected if the magneto- 
spheric boundary were interchange unstable. The Rayleigh- Tailor instability of the boundary can 
occur in bright long-period X-ray pulsars. This instability would be, however, suppressed by the 
magnetic pressure gradient if the X-ray luminosi t y of the neutron star undergoing spherical accre- 
tion is Lx < 3 x 10 36 ergs -1 ( Arons <fe Leal Il976l : Eisner &: Lamblll977l ). and by the magnetic field 
shear in the magnetopause if the spin pe riod of the neutron star is smaller than a few hundred sec- 
onds (see e.g. llkhsanov Pustilniklll996l . and references therein). The Kelvin-Hclmholtz instability 
can be effective in short-period X- ray pulsars in which the mag netospheric radius of the neutron 
star is close to its corotation radius (jBurnard. Lea &: Aronslll983l ). But this instability might not be 
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effective if the relative velocity between the magnetosphere and the accretion flow e xceeds the speed 



of sound or /and if the material at the bou ndary is magnetized (for discussion see lAnzer &: Borner 



198G . Il983l : iMalagoli. Bodo Sz Rosnerlll996l ). This indicates that the interchange instabilities of the 
boundary can be at work only in a limited number of pulsars, while Bohm diffusion takes place 
in any of the considered objects. Having this in mind, we focus our consideration on the case of 
Bohm diffusion assuming that the interchange instabilities of the boundary are suppressed. 

The stationary accretion picture implies that the rate of plasma diffusion into the pulsar 
field at the magnetospheric boundary is equal to the mass capture rate by the neutron star from its 
environment at the Bondi radius and to the mass accretion rate onto the surface of the neutron star. 
Let us assume that this condition is satisfied at the distance r mca , i.e. 97ti n (r mca ) = LxRns/GM ns . 
Combining Eqs. (l28l) . (|29j) and (i30j) . and solving the above equation for r mca , one gets 

r mca c 8 x 10 7 a^ 13 4 13 T 6 - 2/13 m 5 / 13 i?- 4 / 13 L3 7 4 / 13 cm (31) 

Here ao.i = a/0.1, Tq = Tn/10 6 K is the plasm a temp erature at the magnetospheric boundary, 



which is normalized according to iMasetti et al.l (|2006l ). and L37 is the X-ray luminosity of the 
pulsar in units of 10 37 ergs -1 . 

The value of r mca under the same conditions is smaller than the value of r a , which represents 
the canonical radius of the magnetosphere of the neutron star undergoing spherical accretion. This 
can be associated with accumulation of material at the inner radius of the slab which occurs if 
SOtin < 9Jt c . As the material accumulates the gas pressure increases and the slab approaches the 
neutron star to a closer distance. Since Wl- m oc r m 13 ^ 4 the diffusion rate of material into the 
magnetosphere also increases and reaches 5DT C as the inner radius of the slab decreases to r mca . 

The spin-down torque applied to the neutron star from the magnetic slab can be evaluated as 
K^f = fc t 9j?u; s r mca . Taking into account that Wl = [i 2 / (2GM m rl lca ) 1/2 one finds 

^(si) = k t [i 2 uj s = k t fi 2 ^ 

'"mca '"cor ) 

Thus, the spin-down torque applied to the neutron star accreting material from the magnetic slab 
is a factor of ~ (?"a/ r mca) 3//2 higher than the maximum possible spin-down torque applied to the 
star accreting material at the same rate from a quasi-spherical flow. 



7. Signs of magnetic accretion in GX 301 2 

The maximum possible spin-down rate of the neutron star in GX 301-2 within the magnetically 
controlled accretion scenario can be evaluated as z^f? 1 = iffj /2irl, that for the pulsar parameters 
is 

4 mca) ~ TxlO^Hzs- 1 & t a" 3713 m~ 14 /i3 (33) 
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T -l rp3/13 r 6/13 R 57/13 p _i / g* \ 

J 45 J 6 ^37 ^6 ^685 I R I 

\ -DCRSF / 



17/13 



Hence, the spin-down rate of the neutron star observed during the spin-down trends can be ex- 
plained in terms of the magnetically controlled accretion scenario provided the surface field of the 
star is close to -Bcrsf an d the efficiency coefficient is k t > 0.14. The spin-down power of the neutron 
star within this scenario is spent in the energy of electric currents and turbulent motions excited 
at the magnetospheric boundary and, possibly, to spin-up the material at the inner radius of the 
magnetic slab. In this case the angular momentum lost by the star can partly be accumulated in 
the material surrounding its magnetosphere. The spin-down torque in this case decreases as the 
angular velocity of the material at the inner radius of the slab approaches the angular velocity of 
the star itself. The angular momentum accumulated in the slab can later be transferred back to 
the star due to accretion process. The star in this case is expected to switch from the spin-up to 
the spin-down phase without any significant change of the accretion flow geometry. 

A detailed study of the angular momentum exchange between the star and the slab is beyond 
the scope of this paper. Here we would like to note only that the angular velocity to which the 
material in the slab could be spun-up within this scenario is limited to the angular velocity of the 
star itself, oo$. The angular momentum transfer from the star to the slab can bring its inner part to 
a solid rotation. The radial size of this part is limited, however, to the corotation radius at which 
uj s = Wk(^cor)- This opens a possibility to consider a situation in which the spin-up of the material 
at the inner radius of the slab results in temporary transition of the accretion picture from the 
non-Keplerian slab to Keplerian disk of the radius r cor . The star in this case will experience a rapid 
spin- up on a timescale of the vis cous time at the c orotation radius. It is interesting that the size of 



the Keplerian disk evaluated bv lKoh et al.l (|1997l ) from modelling the rapid spin-up events is very 



close to the corotation radius of the pulsar. 



8. Discussion 

One of the basic assumptions adopted in the magnetically controlled accretion scenario is that 
the magnetic field energy density in the material captured by the neutron star at the Bondi radius 
is comparable to its thermal energy density. The latter can be evaluated by taking into account 
that poo = 9Jl/irRQV ie \. This yields 

E* * 0.02 erg cm" 3 m~ 2 c\ M 17 ( g^g^ ) ' • (34) 

The magnetic energy density in the captured wind depends on the field strength of the massive 
companion of the neutron star. The dipole approximation (B mc oc a~ 3 ) to the magnetic field of 
the massive stars remains valid up to a distance a^, at which the dynamical pressure of the wind 
ejecting by this star reaches the magnetic tension of the stellar dipole field (here a is the distance 
from the massive companion). The magnetic field in the material propagating beyond ak decreases 
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as B oc a 2 ( Walder. Folini Mevnetll201ll ). Therefore, the magnetic energy density in the wind 
at the binary separation clq > is £ m (ao) = £m0 = mL/ {^ 7Ta k a o}' which implies 



0.33 erg cm a 



-3 „-4 
13 



10 39 Gcm 3 



100 E(7 



-2 



(35) 



Here /i ms is the dipole magnetic moment of the massive star and 013 is the binary separation in 
units of 10 13 cm. It, therefore, appears, that the magnetic accretion scenario can be expected in 
GX 301-2 if the surface field of the massive component is a few hundred Gauss. 



Recent spectropolarimetric observations (see, e.g. lHubrig et al.ll2006l ; IOksala et alJl2010i ; lMartins et al 



2010i . and references therein) have shown a relatively strong magnetization of O/B-type stars to be 
not unusual. The strength of the large-scale field at the surface of these objects has been measured 
in the range ~ 500 — 5000 G, and in some cases beyond 10 kG. In this light, the assumption about 
relatively strong magnetization of the massive star in GX 301-2 seems to be rather reliable. It 
can be tested by spectropolarimetric observations of this system, which, therefore, appears to be 
of great importance for justification of the magnetically controlled accretion scenario used in our 
paper. 



9. Conclusion 

We have shown that the rapid spin-down of the neutron star observed during the spin-down 
trends of GX 301-2 can be explained in terms of magnetically controlled accretion scenario provided 
the surface field of the star is close to the value, -BcRSFj derived from observations of the cyclotron 
resonance scattering feature in the X-ray spectrum of this pulsar. The same strength of the magnetic 
field has been evaluated from modeling of the rapid spin-up events in terms of the Keplerian 
disk accretion. We conclude that the pro blems in modeling of the spin evolution of GX301- 



2 encountered by iDoroshenko et al.1 (|2010i ) indicate an oversimplification in the currently used 
accretion scenarios rather then an extremely high magnetization of the neutron star. We have shown 
that these problems can be avoided by incorporation of the magnetic filed of the accretion flow into 
the quasi-spherical accretion scenario provided the relative velocity of the neutron star satisfies the 
condition v CT < v rc \ < v mca - The accretion flow in this case is decelerated by its own magnetic 
field at the Shvartsman radius, which exceeds the magnetospheric radius of the neutron star, and 
accumulates in the non-Keplerian magnetic slab confined by the magnetic field of the flow itself. 
The plasma approaches the star on the timescale of field dissipation in the magnetic slab, which 
significantly exceeds the dynamical (free-fall) timescale. The process of plasma penetration into the 
stellar magnetosphere is governed by the magnetic reconnection and drift-dissipative instabilities 
and occurs at a rate of the Bohm diffusion. The magnetospheric radius within this scenario, r mca , 
is smaller than the canonical value, r a . The spin-down torque applied to the neutron star from 
the slab is -fT s d = kt fj?/ (r mca r cor ) 3//2 , where 0.1 < k± < 1. An exchange of the angular momentum 
between the star and the slab can lead to variations of the spin period of the pulsar. 
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Fig. 1. — Pulse frequency and 12-50 keV pulsed flux observed with the Fermi Gamma-ray Burst 
Monitor detectors. GMB Pulsar Project |http:/ /gammaray.nsstc.nasa.gov/gbm/science/pulsars/ 



